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Abstract 

Transition states or quantum states of zero energy appear at the boundary between the discrete 
part of the spectrum of negative energies and the continuum part of positive energy states. As 
such, transition states can be regarded as a hmiting case of a bound state with vanishing binding 
energy emerging for a particular set of critical potential parameters. In this work we study the 
properties of these critical parameters for short range central potentials. To this end we develop 
two exact methods and also utilize the first and second order WKB approximations. Using these 
methods we have calculated the critical parameters for several widely used central potentials. 

The general analytic expressions for the asymptotic representations of the critical parameters 
were derived for cases where either the orbital quantum number / or the number n of bound states 
approaches infinity. 

The above mathematical models enable us to answer the following physical (quantum mechanical) 
questions: i) what is the number of bound states for a given central potential and given orbital 
quantum number I; ii) what is the maximum value of / which can provide a bound state for the 
given central potential; iii) what is the order of energy levels for the given form of the central 
potential. 

It is revealed that the ordering of energy levels depends on the potential singularity at the origin. 
PACS numbers: 03.65.Ge, 03.65. Sq, 03.65.Ta 
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I. INTRODUCTION 



It is generally agreed that estimating the number of bound states of the Schrodinger 
equation is a problem of great practical importance. A substantial effort was devoted for 
evaluating the upper and lower limits on the number of bound states for a given central 
potential. V. Bargmann [l| and J. Schwinger {2] seem to be the first who tackled this 
problem. Since then many authors have contributed to the study of this problem. Among 



them we would like to refer to 
of F. Calogero and F. Brau 



le works js-T] and to emphasize particularly the contribution 



12|. 



The purpose of the current contribution is to present a systematic approach to the in- 
vestigation of quantum states with zero energy. These transition states appear between the 
discrete part of the spectrum of negative energies and the continuum part of positive ener- 
gies. We study transition states by solving the proper Schrodinger equation for short range 
central potentials possessing specific sets of critical parameters. To this end, we develop two 
exact methods for solving the zero energy Schrodinger equation, and for obtaining the values 
of the associated critical parameters. We apply these methods to find the critical param- 
eters of an improtant class of central potentials including among others the Gaussian and 
the Yukawa interactions. A few examples of potentials admitting the analytic solutions for 
the transition states are exhibitted and the associated critical parameters are presented in 
analytic form. We provide numerical results in the form of tables of the critical parameters. 
These results enable one to obtain the exact number of bound states for these potentials 
without any additional computations. 

Analyzing our results we have observed the following universal properties of the solu- 
tions of the Schrodinger equation with central potentials: (a) For a given orbital angular 
momentum quantum number / the n*'^ critical value of the universal parameter Pn,i of any 
potential behaves as for large n. (b) For a given number n of bound states the critical 
parameter Pn,i grows as l"^ with increasing I. Both results can be explained with the help of 
the WKB approximation, (c) The ordering of the energy levels with various {n, 1} depends 
on the potential singularity at the origin. 

The paper is organized as follows. In Section |TT] we briefly discuss the Schrodinger equation 
and the properties of short range potentials. The asymptotic behavior of the solution at 
the origin and at infinity is discussed in sections IllipVt respectively. In section |V] we 
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introduce a phase- kind equation for the transition states, and in section |VT] we describe the 
WKB approximations for calculation of the critical parameters. The analytic asymptotic 
expressions for the critical parameters are presented in section IVIII Numerical results and 
conclusions then follow in sections IVIIII and IIXI 



II. SHORT RANGE POTENTIALS 



Let us start this section with a quotation from section §18 in the Landau & Lifshitz 
textbook on quantum mechanics |13|: "// the field diminishes as —\jr^ at infinity, with 
s > 2, then there are not levels of arbitrarily small negative energy. The discrete spectrum 
terminates at a level with a non-zero absolute value, so that the total number of levels is 
finite". 

We notice, that the statement expounded in the first sentence of the above quotation is 
doubtful. It contradicts our practical experience and also section §133 of the same book 



13| . The conclusion made in the second sentence proves to be true. 

Therefore, we shall consider central potentials V{r) satisfying the corresponding boundary 
condition at infinity: 

lim rV(r) = 0. (1) 

r— >oo 

Near the origin, the boundary condition for an attractive potential can be written in a similar 
form 

limrV(r)=0, (2) 

r-±Q_ 



to avoid fall of a particle to the center [13| (§§18,35). 

We shall consider the class of potentials which can be presented in the form: 

V{T,r,) = -^f(-\ (^7>0, ro>0), (3) 



where g is the coupling constant which determines the strength of the interaction. Note that 
according to condition ([2]) the power s must satisfy the inequality 

g = s+p<2, (4) 

where p corresponds to the leading term in series expansion of function /(r) near the origin 

/(r) ~ ^. (5) 



It should be realized that the definition ([3]) describes a wide class of potentials such as, e.g., 
square well, exponential, Hulthen, Gaussian, Yukawa, Woods-Saxon and many others. 
The radial part of the Schrodinger equation for a single particle moving in a central 



potential field takes the form 



13|(§32): 



X 



0, 



(6) 



where m is the reduced mass, and x = xK'") is the reduced radial part of the wave function 
for a stationary state with angular momentum / and energy E. 

Changing the potential parameters the appearance of a new bound state is accompanied 
with a new solution of Eq.([6]) for E = {). As we limit our discussion to potentials of the form 
([3]), we rewrite Eq.([6]) for zero energy in the form 



2m g 



/ 



+ 



(7) 



The scale transformation x = r/rg leads to the equation: 



dx"^ 



2m gr^ 



2-s 



X° 



-f{x) + 



1(1+1)' 



X^ 



X{x) 



(8) 



It is seen that the potentials V{r,ro) and rg ^V{r,l) have equivalent solutions of Eq.(l6]). 
Thus, we are interested in solving the equation 



X"{r) = [/(r)x(r). 



where 



U{r) = -f3v{r) + 



vir) 



f{r)r~ 



(9) 

(10) 
(11) 



Now our choice of the form ([3]) for central potentials becomes clear. The solutions of Eq. ([9]) 
for a given angular momentum quantum number / depend effectively, (see, e.g.,[l^) only on 
one parameter 

(12) 



_ ^mg . 



Eq. ([9]) is a differential equation of a second order. In order to solve it both analytically and 
numerically one needs to know the behavior of its solution near the origin and at infinity. 
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III. THE SOLUTION NEAR THE ORIGIN 



At first, let us consider the solution of Eq.Q near the origin. Expanding the potential 
into a power series and keeping only the leading term we obtain 



X"ir) 



(13) 



where 

and q is defined by Eq. 
condition 

has the form 



A = 137,, (14) 
The particular solution of Eq. flT^ . satisfying the boundary 



x(o) = 0, 



(15) 



(16) 



where Ju{z) is the Bessel function of the first kind, and A is an arbitrary constant. Keeping 
the first two terms in the series expansion of the function f ll6p . one obtains 



X(r) ~ A r'+^ [(2 - g)(2/ + 3 - g) - A r^-^l . 

r— >0 



(17) 



In general, for integer and half- integer q the solution of Eq.(|9]) satisfying the boundary 
condition ( !T5|) can be presented by the following infinite series 



x(r) = A r 



i+i 



i=l 



(18) 



For integer q all of the 6-coefficients should be zero. It follows from Eq. (fT7|) that 

A 



Cl 



(19) 



(2-g)(2/ + 3-g) 

Substituting the representation ( fTSll into Eq.Q, and then equating the expansion coefficients 
of the same powers of r for the left-hand (Ihs) and right-hand sides (rhs) of Eq.([n]), one can 
calculate any finite number of the subsequent coefficients q and 6j with i > 1. 



IV. THE ASYMPTOTIC BEHAVIOR OF TRANSITION STATES 



For an eigenfunction \1/ belonging to the discrete part of the spectrum, the integral 
J|\I^p(iy, taken over all space, is finite. This certainly means that |\l'p decreases quite 



rapidly, becoming zero at infinity. The motion of tlie system is limited to a finite range, and 
it is said to be in a bound state. For wave functions belonging to the continuous part of the 
spectrum, the integral J\'i/\'^dV diverges due to the fact that |\E'p does not become zero at 
infinity (or becomes zero insufficiently rapidly). 

On the other hand [13] (§ 18), the spectrum of negative eigenvalues of the energy is 
discrete, i.e. all states with E < in the field which vanishes at infinity are bound states. 
The positive eigenvalues -E > 0, on the other hand, form a continuous spectrum. 

In other words, for the bound states the eigenvalues of the energy E < 0, and the 
eigenfunctions must satisfy the boundary condition 

lim 1^1^ = 0. (20) 

In contrary, for a free state that belongs to the continuous spectrum the energy eigenvalues 
E > and |\E'p does not become zero at infinity (or becomes zero insufficiently rapidly). 

For the transition states with E = the asymptotic behavior (r — )■ oo) of the eigen- 
functions remains unclear. 

The important step is to realize that the boundary condition fl20|) must be valid 
for the transition states {E = 0) as well. Thus, for these states achieves zero at 
infinity. However, we note that it may tend to zero too slowly to ensure the convergence of 
the integral f\'^\'^dV. 

In the following we shall rely on the boundary condition fl2U]) for the transition states. 

For / > the asymptotic boundary condition ([T]) enables us to neglect the potential 
V{r) = —I3v{r) in Eqs.(p|)- (fTT|) at large enough r. The general solution of the resulting 
equation has a form 

Xi{r) = Cy^' + C2r-\ (21) 

where Ci and C2 are arbitrary constants. As \I' ~ R{r) = Xii.''^)/^) o^ie should put Ci = 
in order to satisfy the asymptotic condition (1201) . Thus, we get 

Xi{r) ~ C^r-K (22) 

Expressing the latter equation in the form r''xi{r) — const, we obtain the following con- 
dition for the first derivative: 
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It is clear that the asymptotic behavior of the solution xi{^) of Eq.(|9]) depends on the 
parameter (3 of the effective potential ffTOj) . Thus, according to Eq. f l23|) . the solution xi of 
Eq.dn]) fulfills the asymptotic condition 



for the critical parameters Here n is a number of zeros of the function for the 

given potential f lTOj) . Hence, by definition, if for a given I the potential V{r) is characterized 
by the parameters meeting f3n+i > P > Pn then the proper number of bound states equals n. 

The asymptotic condition fl24|) was derived assuming that the orbital quantum number 
/ > 0. However, it is easy to show that Eq. (l2^ preserves its validity also for Z = 0. A typical 
graph of the function is presented in Fig. [T] 

The straightforward solution of the second order differential equation ([9]) with the bound- 
ary conditions ( !T7|) and fl2^ presents our first method for calculating the critical param- 
eters Pn of a given attractive potential ([3]) satisfying the boundary conditions ^ and ([2]). 
This method is especially effective and accurate for small values of /. For a few potentials, 
such as the exponential, the Hulthen and the Woods-Saxon, one can derive analytical ex- 
pressions for the critical parameters using this method (see Appendix). However, this is 
possible only for ^-states (/ = 0), when Eq.([9]) with the potentials mentioned above has a 
general analytical solution. We are familiar with only one form of central potential which 
admits an analytical solution of Eq.(|9]) for / > 0. It is a cut-off potential (described in the 
Appendix) for which the finite square well potential presents its particular case. 

Let us add one important comment. We refer to Eq.( l24l) as the asymptotic behavior 
condition. However, the solution ( 12T|) . and therefore - condition ( 12^ . correspond to the 
assumption that the potential V{r) is negligible in comparison with the centrifugal term 
/(/ + l)/r^. Therefore, the condition ( 12^ is applicable at a distance r when the condition 



Fi{Pn)= lim -x,(r) + x;(r) =0 



(24) 



|y(r)|« 



/(/ + !) 



(25) 



is satisfied. 



V. EQUATION OF THE PHASE KIND 



It was mentioned in the preceding section that the straightforward method for calculat- 
ing the critical parameters of central potentials loses its accuracy with increasing angular 
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momentum quantum number /. In this section we propose another method for calculating 
these parameters. This method is based on the logarithmic derivative y{r) = x'{^)/x{^) of 
the reduced radial wave function introduced earlier. The final equations are close to, but 



10 



12|. 



differ from the so called phase equations presented in 
Let us start with the trivial identity 

^] = ^-(^] , (26) 
Xj X \xj 

and transform the radial Schrodinger equation into a Riccati type equation for the 
corresponding logarithmic derivative y{r) = yi{r): 

y\r)+y\r) = U{r). (27) 

The asymptotic behavior of the logarithmic derivative for / > 

yi{r) - -- (/>0) (28) 

r— s>oo r 



follows from the asymptotic representation ( l22ll . 

To deduce the asymptotic behavior of the logarithmic derivative for the transition S- 
states (/ = 0) let us start with the fact that for this case the rhs of Eq. fl27|) equals V{r). Let 
us then consider central potentials with the asymptotic behavior 

V{r) -A > 0) (29) 

where /i > 2 according to the boundary condition ([1]). The general solution of the proper 
Schrodinger equation ([9]) has a form: 

X{r) = [CsJ. (2iy^r^^ + C^J-, (2u^r^^'^ , (30) 

with 

z/^^<0 (/i>2). (31) 

It is seen that for /i > 2 the argument of the Bessel function goes to zero as r — )■ oo. Thus, 
using series expansion for the Bessel functions, it is easy to show that one should put C2 = 
in order to satisfy the boundary condition (!20l) . Taking then the logarithmic derivative for 
the resulting x(r), and once more using a series expansion for the Bessel functions, one 
obtains: 

Voir) ^ -^r'-". (32) 

r->-oo /i — i 
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It is clear that the asymptotic solution (132|) of the Riccati equation (127|) for / = satisfies 
the following inequahty: 

yl{r) « \y',{r)\. (33) 

r— ^-oo 

The asymptotic solution of the Schrodinger equation ^ with the exponential potential 
is presented in the Appendix (see Eq.( ]A20l) ) for / = 0. The corresponding logarithmic 
derivative 

yl'^^^r) ^ /3e-^ (34) 

r— ^-oo 

satisfies the inequality (!33|) . It is easy to show that the asymptotic representation (r — oo) 
of the Hulthen and the Woods-Saxon potentials reduces to the exponential forms —Pe~^ 
and —f3xQ^e~^, respectively. Hence, the asymptotic solution of the corresponding Eq. (l27|) 
can be presented by the rhs of Eq.f lMl) for the Hulthen potential, and by (3xQ^e~^' for the 
Woods-Saxon potential. The latter logarithmic derivatives certainly obey the inequality f p3|) 
as well. 

It is reasonable to suggest that inequality fl33|) is valid for all of the short range potentials 
(may be excluding only the cut-off potentials). In this case, one can neglect the square of 
the logarithmic derivative in the Ihs of Eq. fl27|) . The solution of the latter equation with 
/ = can be obtained then in the explicit form 

POO 

yo{r) ~ /3 / v{r)dr. (35) 

r->oo 

For potentials with the asymptotic behavior (12^ . formula f l35p gives the asymptotic repre- 
sentation (!32|) . For the exponential potential, the rhs of Eq. (l35|) leads to 
For the Yukawa and Gaussian potentials Eq. (l35|l yields: 

2 

Vo'-'ir) ^ y^^-ir) ^ f^. (36) 

For deriving the latter expressions we used the leading terms of the asymptotic expansions 
of the incomplete gamma function r(0, r) and the complementary error function erfc(r) 
obtained as the results of integration in Eq. (!35l) . 

It is easy to check that the asymptotic logarithmic derivatives ( 136|) satisfy the inequality 

(El. 

The substitution 

yiij-) = /C(r) cot //(r) (37) 
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enables us to transform Eg. (1271) into the following equation for the phase function ri{r): 

v'ii^) = Trr^rrr sm2r]{r) + /C(r) cos^ r]{r) - — — sin^ r]{r). (38) 
2/C(r) /C(r) 

The stabilizing function /C(r) can be chosen in a sufficiently arbitrary manner. The simplest 
choice is /C(r) = 1. In this case, the boundary conditions fl25]) and ([35D define the following 
asymptotic condition for T]{r): 

r]{oo) = -7i/2 + n7r. (n = l,2, ...) (39) 

Any function /C(r) that preserves the limit, linir^oo [y{f)/^{^)] = provides the asymptotic 
behavior (1391) . It was established (at least numerically) that the more precise results are 
provided by stabilizing functions of the form: 

JCi{r) = (40) 

where Ui{r) governs the behavior of f/(r) both near the origin and at infinity, that is 

Pv(r), I = 

Ui{r)={ " . (41) 

/(/+l)r-2, />0 



Use of the stabilizing functions ( l40|) - fl4T|) enables us to replace Eq. fl38|l with two simpler 
equations: 

sin 2r^{r) + v^Z^W, = 0) (42) 



r]'(r) = - + 1) cos2?7(r) sin2r7(r 

r 2 



:sin2r/(r). (/ > 0) (43) 



Substituting expressions ( l40l) -( l4T]) for / > into definition ( 1371) with the asymptotic form 



( l28l) . one obtains: cot 77(00) = —^Jl/{l + l). Thus, the critical parameters /3„ must provide 
the following asymptotic behavior for the function 17 (r): 

Fi{(5n) = r?(oo)|^=^„ = 5, - ^ + nvr, (n = 1, 2, ...) (44) 

with 



61 = arctan ^ [ i ' 

Eq.dH]) was derived for / > 0. It describes the asymptotic behavior of the solutions of 



Eq. (j43l) . However, it is clear that the stabihzing function ICo{r) = ^/Uo{r) preserves the 
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correctness of the asymptotic formula (139|) for the logarithmic derivatives yoir) with the 
asymptotic behavior defined by Eq. fl35|) . Therefore, condition fH^ - fj45l) with / = can be 
used also for the asymptotic approximation of solutions of Eq. ( 



It is seen from Eq.f H2]) that the stabilizing function /Co(^) = \/Uo(jj can be applicable 
only in case of its nodeless character, otherwise, the simplest choice is ICo{r) = 1. 

The technique described above in this section is sufficient for presentation of the sec- 
ond method for calculating the critical parameters of central potentials of the from (|3]). 
However, we would like to make some additional remarks that can be useful. 

Eqs. fH2]) . f H5]) . along with the boundary condition fHl|) - fH5]) . provide a stable and accurate 
solution to the problem of critical potentials for both small and large values of /. 

A typical graph of function t = Fi{P) has a staircase form. It is presented in FigEl The 
abscissas of the points of the staircase function intersections with lines t = 6i — ti /2 + im 
give the desired critical parameters /3 = /3„. 

As an additional useful information it can be shown that 

Fi{0) = 6i. (46) 

This is because of the following. Setting (3 = 0, the second term disappears from the rhs of 
Eq. fj42l) . The analytic solution to the resultant equation has a form: 

r]{r) =^ arctan (^C^/Uo{r)^ . (47) 

Condition ([T]) thus provides ^^(oo) = for /3 = and I = according to Eq. fHTl) . 

Putting /3 = in Eq. (H3l) . the latter loses the term with f/o(r). The resultant equation 
has analytic solution of the form: 

(2/ + l)tanh(^^^^^lnr + C) -1 



r](r) = arctan , . (48) 

/3^o 2v//(/ + l) 

For arbitrary finite real C, tank presented in Eq. (l48l) approaches 1 as r — )■ oo. The resultant 
expression thus reduces to 6i. 



VI. THE FIRST AND SECOND ORDER WKB APPROXIMATIONS 

In this section we apply the first-order and the second-order WKB approaches to cal- 
culating the critical parameters of central potentials. Unlike the exact methods presented 
earlier these methods are certainly approximate but they are also much simpler. 
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Specific modification of the first-order WKB approach presented below cannot provide 
results of high precision. However, its accuracy grows rapidly with increasing the number n 
of bound states. It is important to note that this method can be applied to transition states 
with any orbital quantum number Z > 0. 

The accuracy of the second-order WKB method grows with increasing /. For example, 
the relative error for I > 10 can be less than 10~^. This enables one to test the results 
obtained by the phase kind method (Sec|V]) for large /. On the other hand, this method can 
be used by itself if a very high accuracy is not needed. The disadvantage of this method is 
its inapplicability for the transition S-states (/ = 0). 

For E = the presence of the centrifugal term in the effective potential U (r) ensures the 
existence of two turning points for attractive potentials of the form ([S]). In the second-order 
approximation [l^ the WKB quantization condition 
can be written as 



17H19| , as applied to our consideration. 



(n 



1,2,...) 



(49) 



where n is the number of bound states as in the previous sections. The term 



^0= / V-U{r)dr (50) 
7r/2 correspond to the first order WKB-approximation (for 



together with the term 5*1 = 
two turning points). The turning points ri and r2 are the roots of the equation 

/(/ + 1 



1,2 



(51) 



where the function t>(r) is defined by Eqs.Q, (fTT]) . The second-order correction takes the 

[ ! 

form 



^2 



lim , 

Ai^+o V 48 



U"ir) 



-dr 



-Xbi\ai\ 



-3/2 



&2|a2| 



-3/2- 



(52) 



where ai, bi, a2 and 62 are the expansion coefficients of the effective potential U{r) in the 
neighborhood of the turning points. That is 



U{ri + 5r) = ai6r + hi{5rY + 



(53) 



where from Eq. (!5T]) it can be seen that 



-/3 



—v{ri) + v'in) 
6 



(54) 
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Notice that the small magnitude of the second-order correction (152|) results from the differ- 
ence of two large terms. Therefore, both these terms must be calculated with high accuracy. 
Nevertheless, we would like to emphasize that this correction increases the accuracy by two 
to three orders of magnitude. 

A modification of the first-order WKB method, where the centrifugal potential is ex- 



cluded from the quasiclassical momentum 28|, can be applied for calculating the critical 



parameters under consideration. According to this approach, the quantization condition for 
the transition states {E = 0) reduces to: 

a/ I3v{r)dr = 7i{n + 7^,5), (55) 



where 



21-1 
4 



q<0 



li,,= { 2/-l + g • (56) 



< g < 2 



2(2 -g) ' 

In accordance to its definition the parameter q is ruled by the behavior of the central 
potential near the origin. The cut-off potential f lAip presents an exclusion. The exact 
analytic solution for this case is presented in the Appendix. 



VII. THE CRITICAL PARAMETER ASYMPTOTICS 

It was shown in the previous sections that in general the critical parameters for central 
potentials can be calculated numerically. For a few special cases, presented in the Appendix, 
one can deduce the analytical results. 

In this section it will be shown that it is possible to derive the analytical expressions for 
the asymptotic form of the critical parameters (3 = /3n,i- In so doing, asymptotic implies a 
situation where either the number of bound states n approaches infinity for a given finite 
/, or the orbital angular momentum quantum number / goes to infinity for a given finite n. 
The first order WKB approach is applied to solve this problem. 

First, let us consider the case of finite n where Z — ?■ 00. 

Numerical calculations demonstrate that the distance between turning points |ri — 
reaches zero as I approaches infinity. This result can be explained and supported by the fol- 
lowing arguments. According to Eqs. (l49l) . (l50l) the standard Bohr-Sommerfeld quantization 
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condition for the transition states reads 

^/-U{r)dr = {n - i)7r. (57) 

This means that the effective potential U{r) must be negative in the range [ri,r2]. Hence, 
according to Eq. ffTOl) the critical parameter /3 must tend to infinity as or faster for / — oo. 
The latter in turn implies that the integrand in the LHS of Eq. (!57|) approaches infinity 
as I — 7- oo, whereas the RHS of Eq. fl57|) remains finite. The above contradiction can be 
eliminated only by setting ri = r2 for the limits of integration, which proves the statement. 

It is clear that the point Rm of minimum of the effective potential U{r) is localized in 
the region [ri,r2], that is 

ri<Rm<r2. (58) 

Therefore, for / approaches infinity Rm tends to the point where the turning points ri and 
r2 merge. Hence, in order to calculate = hm;_s.oo Rm it is enough to solve the following 
set of equations: 

f/(rj = ^ ^ 

(59) 

U'irm) = 

Substituting the explicit form (ITOl) into (1591) and eliminating /3 and /, one obtains the fol- 
lowing simple equation for rm'- 

2virm) + Tmv'irm) = 0. (60) 
Now, any of the two equations fl59|) gives the required asymptotic expression: 

Pn,l^ f/(/ + l), (61) 

with 

= 2 ^ V (62) 
The second derivatives di = d'^[3n,i/ dt^ of the asymptotic critical parameters fl6T|) are 
presented in Table |I]along with the values of ; = — '2Pn,i + Pn,i-i which approximate 
the general second derivative d"^ jUn^i/ dt^ numerically. It is seen from this Table that the values 
of di are very close to the values of Ai 19 for the lowest (nodeless) energy states {n = 1) and 
/ = 19. It is worth noting that even though the Yukawa and Gaussian potentials are very 
different, their asymptotic behaviors are coincident {di = 2e). 



14 



TABLE I: The asymptotic representations of critical parameters /3„ / — ttK^ + 1) ^^'^ 



/3n,/ - + 7/,a)^- The numerical results for A„ ; = /3n,i+i - 2f3n,i + f3n,i-i and 

^n,z = /3ra+i,« — '2f3n,i + Pn-i,i are presented for comparison. 



potential 


exponential 


Hulthen 


Yukawa 


Gaussian 


WS(xo = 1) 


WS(xo = 0.001) 


v{r) 


e~' 


(e^- - 


e~°'' jr 


e 


(l + xoe^-^ 


(l + xoe'^)-i 


fm 


2 


1.59362 


1 


1 


2.21772 


6.17241 


di 


^/3.69453/ 


3.08828 


2e/5.43656/ 


2e/5.43656/ 


4.14224 


0.077658 


Al,l9 


3.69449 


3.08826 


5.43654 


5.43639 


4.14220 


0.077664 


^5,19 


3.69674 


3.08845 


5.43698 


5.43895 


4.14471 


0.077846 


/o°° \/v{r)dr 


2 


vr 












4/4.93480/ 


2 


7r/3.14159/ 


47r/12.5664/ 


6.35257 


0.286909 


Al9,0 


4.93480 


2 


3.14127 


12.5694 


6.35257 


0.286908 


Al9,5 


4.88254 


1.98432 


3.11141 


12.5151 


6.29585 


0.286015 


^-99,5 


4.930 


1.998 


3.138 


12.561 


6.347 


0.2868 



Notice that a set of equations ( 1591) for calculating cannot be applied to the cut-off 
potential ( 1A1I) . The correct result ( 1A24P can be obtained by setting = 1, that is by 
equating the merging point and the matching point r = 1. 

Now let us consider the case of n — >■ oo (/ is finite). To this end, one can successfully 



28 



presented in the previous section. According to the 



employ the modified WKB method 
authors of 28|, their method is "exact in the asymptotic limit — )■ oo", where Ur = n — 1. 
Our numerical results confirm this assertion. Thus, using directly the quantization condition 
( l55ll . one obtains: 

/3n,/ ^ -^{n + 7i,gY, (63) 
where ^ 

dn = 2( ) , (64) 

and the parameter 7^^^ is defined in Eq.(l56l). 

For the FSW-like potential Eq.i lMj) gives dn = 27r^ (^)^) which certainly coincides with 
the analytic solution presented in the 1— section of the Appendix. In general, it is easy to 
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check that the asymptotic expressions for the analytic solutions presented in the Appendix 
coincide with the results of this section. 
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VIII. NUMERICAL RESULTS 



The second derivatives dn = (P(3n,i/dn'^ of the asymptotic critical parameters ( 163|) are 
presented in Table |T] along with values of An^i = (3n+i,i — 2f3n,i + Pn-i,i which approximate the 
general second derivative d'^Pn,i/dn'^ numerically. It is seen that for the numerical second 
derivative A^^; the speed ui of a convergence to the asymptotic value di {I oo) depends on 
the number of bound states n. The fastest convergence corresponds to the smallest n = 1. 
A similar situation is observed for the speed m„ of the convergence of An,i to dn {n oo). 
In this case, the fastest convergence corresponds to the smallest / = 0. Moreover, Table |T] 
shows that ui far exceeds m„. 

Using the first and second methods described in Sections IIVI and N ] respectively, we have 
computed the critical parameters for a few widely used central potentials included in the 
nonrelativistic Schrodinger-equation. The results for the exponential, Hulthen, Yukawa and 
Gaussian potentials are presented in Tables HTHVt respectively. The critical parameters for 
the Woods-Saxon potential with Xq = 1 and Xq = 0.001 are exhibited in Tables |VT] and |VTT1 
respectively. The first value of parameter xq presents the minimal value of parameter R = 0. 
By contrast, the second value of parameter xq corresponds to the case of a large value of 
R/ro corresponding, e.g., to the optical-model calculations |20 |. 

The critical parameters for the exponential, Hulthen and Woods-Saxon potentials appear- 
ing in the Schrodinger equation with the orbital angular momentum quantum number / = 
(^-states) can be calculated analytically. For the exponential and Woods-Saxon potentials 
the requested solutions (/ = 0) are proportional to the squares of zeros of the corresponding 
special functions (see Appendix). It is worth noting that the critical parameters for the 
Hulthen potential (/ = 0) have an especially simple form. They are equal to where the 
number of bound states equals n. All these parameters are displayed in Tables HIl IHIt I VI I 
and IVIII for convenient comparing with the cases of Z > 0. Notice that the relative difference 
between the results obtained by the first method and the analytical ones is less than 10~^^. 

The critical parameters for the cut-off potential of the form ( lAll) are not numerically 
presented here, because there is no problem to provide the proper calculations according to 
Eq. (]A7|) for any orbital quantum number / and parameter s. However, we have performed 
the corresponding computations in order to test both the first and second methods. The 
relative difference was less than 10~^^ for the second method (Sec|V]) for / < 20. The first 
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method provides the same accuracy only for small values of / < 3, whereas for large / this 
accuracy can be provided only for I + n < 20. One should notice that all computations 
were performed by means of the simplest Mathematica-7 codes using the standard (default) 
working precision. It is possible, of course, to enhance the calculation accuracy using, e.g., 
the better working precision. 

Due to lack of space we have restricted the results in our Tables to eight significant figures 
and values of Z + n < 16. 

The Yukawa potential is the only one for which we have revealed some earlier results 



on critical parameters [2l|. They are presented there as the critical screening length for 
the one-electron eigenstates which were obtained in frame of standard energy calculations. 
Those results are limited by 5 significant figures and l + n < 9, and coincide practically with 
those exhibited here in Table IVII 

For the Gaussian potential we have found only the results of the binding energy cal- 
culations (see, e.g., Ref. 22I25I). These energies were computed for / + n < 8 and were 
completely consistent with the critical parameters presented here in Table IVIII 



IX. CONCLUSIONS 



In conclusion, we would like to emphasize that the critical parameters are not of some 
particular character. 

First, they present some universal characteristics of central potentials. 

Second, using the tables of these parameters, one can answer the following questions: 

1) What is the number n of bound states for the given central potential and given orbital 
quantum number /. 

2) What is the maximum value of / which can provide a bound state for the given central 
potential, or vice versa, what is the minimum critical parameter which can provide a bound 
state with a given I for the given central potential. 

3) What is the mutual arrangement (order) of the energy levels En^i (characterized by 
the quantum numbers n and /) for the given form of central potential. 

It is clear that the binding energy En,i rises as the number n of bound states or the 
orbital angular momentum quantum number / increases. Tables HTllVIII show that the critical 
parameters Pn,i exhibit the same properties in respect to the numbers {n, I}. It is important 
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to realize that for any two sets {rii, li} and {n2, h} it follows from inequality 

/3ni,Zi > /3n2,l2 (65) 

for the given central potential that 

Eni,h > En2,l2- (66) 

Therefore, from the presented Tables we can deduce the following important properties of 
the discrete energy spectrum of the considered central potentials: 

En,i > En+i,i-i for the Hulthen and Yukawa potentials, 

En^i < En+i^i-i for the exponential, Gaussian and Woods — Saxon potentials. 



This probably relates to the fact that the Hulthen and the Yukawa potentials are singular 
at the origin. 

It was established that the leading terms of the asymptotic expansions of the critical 
parameters f3n,i have the following forms: 

{aif I — )■ oo 

(68) 
a„n^ n — > oo 

where the general analytic expressions for the factors ai = di/2 and a„ = dn/2 are presented 
in section IVIIi The first of the relationships ( l68l) is valid for a finite number n of bound 
states, whereas the second one is valid for a finite orbital angular momentum quantum 
number /. 
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Appendix A 



1. Finite square well like potential 

Let us examine potential of the form: 

V{r) = 



9_ 
0, 



r < ro 

r > Tq 



{s<2) 



(Al) 



For s = one obtains the potential which is widely known as the finite square well. Eq. 
with the potential (lAip reduces to the form: 



X'{.r) 



4^(1 _ ^ !(i±i) 



(A2) 



where H{x) is the Heaviside step function. Parameter /3 is defined by Eq.f ll2l) . Particular 
solution satisfying the asymptotic condition ( l20l) and vanishing at the origin has a form: 



X{r) 



V 



-I 



with 



Br 



2-s 



r < 1 
r > 1 



(A3) 



2/ + 1 



a 



(A4) 



2 ' 2-s 
Here Ja{z) is the Bessel function of the first kind, whereas A and B are arbitrary constants. 
Matching the logarithmic derivatives of solutions flA3|) at the point r = 1, one obtains the 
following equation for /3: 



VP 



V 



lJa+1 



-/(2/ + l)J„ 



(A5) 



Using the properties of the Bessel functions 
simplest equation: 



J 21-1+S ( — ^ = 
\2-s) 



15| . one can reduce Eq. flASP to the following 

(A6) 



The explicit solution of Eq. (1A6P has a form: 



/3. 



2-s 



J21^, n 



(n = l,2,...) 



(A7) 



where presents the n positive zero of the Bessel function J^i{z). 
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From expansion (9.5.12) [15] for large zeros one has: 



V \\ 4Z/2-1 

iun — [ n -\ TT 7 



A{Av'^ - l)(28z/2 -31) 



2 4 



^+o-7rr 192 



z/ 1 



2 4 



(A8) 



Putting 



2 



= (A9) 

and using Eq. (IA7l) . one obtains the following asymptotic expression for the critical param- 
eters corresponding to potential (]A1|) : 

^ f^VvrV. (AlO) 



This yields, in particular: 



dn? 



/3„(s = 0) ~ 27r^ 



Pn{s = l) ^ vrV2. 



(All) 



From expansion for zeros of the Bessel functions of large order 261] (see, also (9.5.14) 15|). 
one has: 



2k: 



(A12) 



The first coefficients ak,n for {A;,n} < 5 can be found in [26|. Thus, using Eq. (IA9p . one 
obtains for large enough I: 



I3n = (3n,l = l- 



1 - S 



1 + ^ ak,n 



k=l 



2l-l + sY 3 
2-s 



(A13) 



It is seen that the leading term of the asymptotic (/ — ?■ oo) expansion of the critical parameter 
is given by 



Pn,l , ^ 



(A14) 



2. Exponential potential 



The exponential potential has a form: 



V{r) = —g exp 



(A15) 
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Eq.(|9]) with the potential f lAlSP reads 



X"{r) 



—(3 exp{—r) + 



X(r). (A16) 



Parameter /3 is defined here by Eq. (!T2|) with s = 0. Eq. (]A16P has an analytical solution only 
for the case of / = 0. Such a particular solution satisfying the boundary condition f|T5|) has 
a form: 

X(r) = A [Jo{q)Yo {qe-^) - Y,{q)Jo (ge-i)] , (A17) 

where q = 2a/^, whereas Jo{z) and Yo{z) are the Bessel functions of the first and second kind, 
respectively. The argument of the Bessel functions in Eg JA17P achieves zero as r — )■ oo. 
Therefore, using series expansion for the Bessel functions [15] 



.2 



Jo(ge-i) - >o(ge-5) ~ -L (A18) 

one should put 

Mq) = 0, (A19) 

in order solution ( lAlTI) satisfies the asymptotic boundary condition ( l20|l . Taking into account 
condition ( 1A19I) for the transition state, and the first of expansions ( 1A18I) one can write for 
the asymptotic behavior of the solution ( ]A17I) : 



X{r) ^ -AYo{2^(3){l-(3e-^). (A20) 



From Eq. (IA19P one obtains for the critical parameters: 

f3n = \jln, (n=l,2,...) (A21) 

where jo,n presents the n^^ positive zero of the Bessel function Jo{z). 
Putting = in expansion flASP for large zeros, one obtains: 

/ 1\ 1 124 , 



This yields 



^ ^L-\)\ (A23) 

n— >oo 4 \ 4 / 



The value of second derivative limn^oo d'^Pn/dn'^ = 7ir^/2 coincides with the corresponding 
value (lAllI) for the FSW-like potential with s = 1. 
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3. Hulthen potential 



For the Hulthen potential 



V{r) 



9 e 



Eq.Q becomes 



1 — e '■0 



f3 + 



1 



For this differential equation one can obtain a general analytic solution of the form 



x{r) = (j){x) = A X " 2-^i(— a, —a; 1 — 2q;; x) + -B 2-^i(a, a; 1 + 2a; x) 



(A24) 
(A25) 

(A26) 



with 



a 



'15, X = e^ (A27) 

only for S'-states (/ = 0). Here 2-^1 c; 2;) is the Gauss hypergeometric function, A and 
B are arbitrary constants. 

Formulas (15.1.20) and (6.1.18) ^ISi] yield: 

r T{a + 1/2) 



2Fi(a, a; 1 + 2a; 1) 



(A28) 



v^r(a + l) ' 

where T{z) denotes Euler's gamma function. The latter representation enables us to obtain 
the vanishing at the origin (r — t- ^ x — i- 1 ) solution in the form: 



X(r) = (p{x) = C 



AY A \ 

— r(l — a)T — ha 2-Fi(— a, —a; 1 — 2a; x) — 
xj \2 J 

x\^ ( \ \ 

-j r(l + a)r(--aj 2-^i(a, a; 1 + 2a; x) 



(A29) 



with arbitrary constant C . For examination of the asymptotic behavior of solution (IA29p . 
one can use formula (15.3.13) [l5|, which yields: 

2r(2a 



-x)° 2-^1 («, a; 1 + 2a; x) 



Inserting the asymptotic representation (]A30I) for a = a and a 
one obtains after some transformations: 



In(-x). (A30) 
= -a into the rhs of Eq. (IA29D , 



On the other hand, series expansion of solution (IA29P near x = 1 (r — )■ 0) yields: 



x(r) ^ — 2Ca^7r^''^ csc(7ra)r. 

) — >-0 



(A31) 



(A32) 
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Putting 



C 



(A33) 



one can get rid of a-dependence for the leading term of the xir) series expansion. Substi- 
tuting expression (1A33P into the asymptotic representation (lA3ip . one finally obtains: 



SmlTTttl 

~ («7r + r). 



Thus, in order to satisfy asymptotic condition fl2U|) . one should put: 



(A34) 



sin(7ra) = 0. 



(A35) 



The roots of Eg. flA35p are the integers, that is q;„ = n. Thus, from definition flA27p one 
obtains that the critical parameters for the Hulthen potential can be determined from the 
simplest relation: 

--n\ (n=l,2,...) (A36) 



^2 

where n is a number of S-bound states. 



4. Woods-Saxon potential 

Finally, let us examine the Woods-Saxon potential 

9 



V{r) 



1 + exp( 



{R > 0) 



which is the most complicated one. For this case Eq.Q takes a form: 



x"{r) 



/3 



+ 



Xir), 



(A37) 



(A38) 



1 + Xoe'" 

where (3 is defined by Eq. f|T2|) with s = 0, whereas Xq = exp{—R/rQ). Eq. flA38P admits 
analytical solution only for the case of Z = 0. Introducing a new variable x = xqc^ and a 
new parameter a = a/^, one obtains a new differential equation 

,2 



a 



X ip"(x) + xip'(x) H vix) = 

1 + X 



(A39) 



for function (p{x) = xi^^)- The vanishing at the origin solution of Eq. flA39P has a form 



(p{x) = C a, XQ)F{a, x) — F{a, xo)F{—a, x)] 



(A40) 
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where 

F{a,x) = x'" 2Fi{ia,ia; 1 + 2ia; -x). (A41) 

Here 2Fi{a,b] c; z) is the Gauss hypergeometric function, C is arbitrary constant. For 
considered case, formula (15.3.13) [l5| yields 

(x)SFi(6,6;1 + 26;-x) ^ ?^^\n{x). (A42) 



Inserting the latter representation into solution (]A40j) and returning to the initial variable 
r, one obtains: 



x(r) ~ 2C(r + lnxo) 



F(— a, a;o)r(2ia;) F{a,xo)T{—2ia) 



2C(r + In XQ)ia 2F1 (-ia, ia; 1; ) . (A43) 



Thus, to satisfy the condition (120|) for the transitional states, one should put 

2F1 ( 1; -— ) = 0. (A44) 



The roots /3„ = of the latter transcendental equation present the desired critical pa- 
rameters for a given xq. It is worth noting that Eg. (1A44I) can be simplified if one uses the 
following relationships between the Gauss hypergeometric functions, the Jacobi functions 
pj)"''''^ and the Legendre functions P,^: 

2Fi(-^a,^a;l;-^) = p/„°-')(2z + l) = (1 + ^)pS(2z + 1) = i?e [Pi„(2^ + 1)] 

(a > 0, z> 0). (A45) 

According to the asymptotic expansion for the Legendre function of imaginary degree (see, 
Eq.(3.2) [27]), one has 

Re [PUcosh t)] ~ i= jyik - l)VMt) (^-^J .h{at) + O {a-""-') , (A46) 

where 

/ , f t\ ao(t) f cosh t - 2 1 

In zero approximation for large enough a, one can put = 0, whence 



Re[Pia{cosht)] ~ t coth(l]Jo{at). (A48) 



V2\l \2 
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Thus, roots of Eq. (lA44p for large enough a = are very close to zeros of the Bessel 
function Jo(ct^), where 

t = arccosh + 1^ = 2arcsinh (^7^^ = 2 In ^ ^ -^ . (A49) 



In terms of zeros jo,™ of the Bessel functions, one can then write down: 

,o,„-2v/^lnfi±^^V (A50) 



Taking into account expansion (IA22p . one obtains finally: 



Pn ^ 7 — ^ , (n--^ . (A51) 



41n^ 
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TABLE II: Critical parameters /3 = 2mgrQh of the exponential potential V{r) = —g exp(— r/ro) 



n\l 







1 



1 
2 
3 
4 
5 
6 
7 
8 
9 

10 
11 
12 
13 
14 
15 
16 



1.4457965 7.0490613 
7.6178156 16.921126 
18.721752 31.525958 
34.760071 50.947660 
55.733076 75.226301 
81.640838 104.38402 
112.48338 138.43438 
148.26072 177.38629 
188.97285 221.24597 
234.61978 270.01793 
285.20151 323.70558 
340.71804 382.31153 
401.16937 445.83785 
466.55550 514.28624 
536.87643 587.65804 473.25639 
612.13217 472.63339 415.36754 



16.312928 
29.879667 
48.076670 
71.002111 
98.713352 
131.24653 
168.62580 
210.86809 
257.98574 
309.98801 
366.88206 
428.67350 
495.36688 
566.96589 



29.258323 
46.518231 
68.345880 
94.837734 
126.05709 
162.04738 
202.83958 
248.45649 
298.91536 
354.22957 
414.40969 
479.46420 
549.40007 



45.892427 66.218077 

66.845206 90.863808 

92.323606 120.00486 

122.41832 153.72542 

157.19351 192.08825 

196.69590 235.14112 

240.96044 282.92106 

290.01388 335.45733 

343.87710 392.77335 

402.56672 454.88817 

466.09616 521.81745 

534.47640 479.64012 



90.236557 117.94852 
118.57542 149.98069 
151.38676 186.46751 
188.74853 227.48127 
230.72116 273.07970 
277.35218 323.30940 
328.67953 378.20847 
384.73391 437.80834 
445.54045 502.13521 
511.11996 488.5090 



149.35431 
185.07997 
225.24597 
269.91954 
319.15566 
372.99967 
431.48931 
494.65627 



15 



14 



474.58293 
415.99092 
361.79612 
13 



476.68376 419.44387 
417.32411 363.76075 
362.41992 312.54336 
311.91912 265.736481 



12 



11 



483.54571 426.4139 
422.43924 368.94348 
365.90248 316.06058 
313.89306 267.72133 
266.36122 223.87344 
223.24818 184.45415 
10 9 
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TABLE III: Critical parameters /? = Imgrgh ^ of the Hulthen potential 



= -g exp(-r/ro)/[l - exp(-r/ro)]. 
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169 


200.99951 


236.28888 


274.80601 


356.85339 


311.46707 


269.16916 


229.95964 


4 


14 


196 


230.24882 


267.81254 


367.49439 


321.39765 


278.38907 


238.46863 


201.63626 


3 


15 


225 


261.49399 


378.38000 


331.57422 


287.85660 


247.22711 


209.68569 


175.23229 


2 


16 


256 


389.50148 


341.98731 


297.56137 


256.22366 


217.97415 


182.81284 


150.73968 


1 




15 


14 


13 


12 


11 


10 


9 
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TABLE IV: Critical parameters /3 = 2mgroh ^ of the Yukawa potential V{r) = —g exp(— r/ro)/r. 



n\l 







1 
2 
3 
4 
5 
6 
7 
8 
9 

10 
11 
12 
13 
14 
15 
16 



1.6798078 9.0819590 

6.4472603 17.744576 

14.342028 29.461426 

25.371660 44.261254 

39.538842 62.160193 

56.84486 83.168247 

77.290455 107.29208 

100.87607 134.53636 

127.60202 164.90453 

157.46853 198.39917 

190.47575 235.02231 

226.62381 274.77554 

265.91281 317.66016 

308.34282 363.67724 

353.91391 412.82768 658.19765 

402.62614 682.98698 599.51605 



21.894984 
34.420414 
49.969576 
68.571467 
90.245270 
115.00434 
142.85836 
173.81459 
207.87862 
245.05485 
285.34681 
328.75740 
375.28899 
424.94360 



40.135552 
56.514114 
75.899394 
98.317925 
123.78892 
152.32675 
183.94242 
218.64457 
256.44012 
297.33466 
341.33284 
388.43851 
438.65494 



63.808976 
84.036777 
107.26037 
133.50366 
162.78498 
195.11871 
230.51632 
268.98711 
310.53874 
355.17758 
402.90900 
453.73757 



92.917164 
116.99234 
144.05569 
174.12875 
207.22876 
243.36968 
282.56299 
324.81824 
370.14347 
418.54550 
470.03018 



127.46092 
155.38248 
186.28656 
220.19272 
257.11705 
297.07295 
340.07168 
386.12280 
435.23449 
487.41381 



167.44064 212.85654 
199.20797 248.46926 
233.95349 287.05673 
271.69505 328.63535 
312.44769 373.21919 
356.22417 420.82039 
403.03541 471.44949 
452.89086 525.11571 
505.79870 



15 



14 



634.12021 
576.31864 
521.48163 
13 



610.77114 

553.83444 
499.87607 
448.88371 
12 



588.17003 
532.08136 
478.98522 
428.86991 
381.72227 
11 



566.34038 
511.08099 
458.82871 
409.57254 
363.30010 
319.99729 
10 



545.31066 
490.8596 
439.43054 
391.01331 
345.59639 
303.16658 
263.70873 
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TABLE V: Critical parameters f3 = 2mgrQh of the Gaussian potential V{r) = —g exp(— r /rg 



n\l 







1 



1 
2 
3 
4 
5 
6 
7 
8 
9 

10 
11 
12 
13 
14 
15 
16 



2.6840047 12.099309 

17.795700 35.089777 

45.573480 70.482856 

85.963400 118.35165 

138.94811 178.72771 

204.51926 251.62783 

282.67201 337.06209 

373.40327 435.03714 

476.71088 545.55758 

592.59331 668.62678 

721.04936 804.24729 

862.07814 952.42108 

1015.6789 1113.1497 

1181.8511 1286.4345 

1360.5941 1472.2765 
1551.9077 



26.901078 
57.675772 
100.71505 
156.12749 
223.97020 
304.27712 
397.07024 
502.36469 
620.17137 
750.49845 
893.35221 
1048.7376 
1216.6587 
1397.1187 



47.107862 
85.627393 
136.32500 
199.31846 
274.67746 
362.44675 
462.65690 
575.32978 
700.48163 
838.12487 
988.26926 
1150.9226 
1326.0912 



72.733004 103.78395 140.26479 182.17789 
118.97711 157.74107 201.92788 251.54239 
177.33962 223.77276 275.63189 332.92115 
247.93811 301.99288 361.48565 426.41761 
330.84800 392.47880 459.56651 532.10794 
426.12056 495.28556 569.93140 650.04986 
533.79230 610.45368 692.62349 780.28811 
653.8902 738.01395 827.67628 922.85815 
786.43489 877.99030 975.11633 1077.7885 
931.44233 1030.4018 1134.9651 
1088.9253 1195.2637 



229.52463 
306.58749 
395.64293 
496.78909 
610.10035 
735.63442 
873.43695 
1023.5447 



1258.8943 



972.06166 



923.18228 826.67993 



750.96802 
713.13988 627.74431 



15 



14 



790.77458 
663.31706 
547.78488 
13 



876.76204 782.42916 693.54140 
832.66745 740.53949 653.85026 572.60004 
700.87689 616.41518 537.38923 463.79866 
581.10160 504.32139 432.97603 367.06499 
473.26149 404.17400 340.52222 282.30589 



12 



11 



10 



9 
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TABLE VI: Critical parameters j3 = 2mgrQh ^ of the Woods-Saxon potential 
V{r) = + exp(r/ro)]. 



n\l 







1 
2 
3 
4 
5 
6 
7 
8 
9 

10 
11 
12 
13 
14 
15 
16 



1.7205730 8.2135317 18.813940 33.542490 52.406813 

9.6742198 20.755782 35.931282 55.231021 78.665326 

23.969284 39.431300 59.073131 82.876856 110.83541 

44.615717 64.332441 88.346797 116.58590 149.01715 

71.614405 95.502170 123.81377 156.42640 193.28014 

104.96555 132.96409 165.51246 202.44471 243.67452 

144.66921 176.73264 213.46842 254.67366 300.23751 

190.72542 226.81735 267.69958 313.13737 362.99747 

243.13418 283.22488 328.218988 377.85407 431.97649 

301.89551 345.96008 395.03646 448.83786 507.19210 

367.00940 415.02662 468.15961 526.09990 588.65845 

438.47587 490.42734 547.59441 609.64918 676.38713 

516.29490 572.16450 633.34572 699.49305 

600.46650 660.23992 725.41749 

690.99068 754.65510 
787.86743 



75.410024 
106.23821 
142.94481 
185.62315 
234.34169 
289.15197 
350.09364 
417.19781 
490.48934 
569.98836 
655.71141 



102.553558 133.83814 

137.9514072 173.80576 

179.2025315 219.60692 

226.3936688 271.32232 

279.5910514 329.01551 

338.8462258 392.73727 

404.2000006 462.52874 

475.6851818 538.42365 

553.3285074 620.44998 

637.1520405 591.78238 



169.26419 
213.80169 
264.15685 
320.40492 
382.60663 
450.81159 
525.06071 
605.38815 



558.22885 



549.06873 483.79929 
540.75901 475.91083 415.20481 
533.21997 468.80008 408.52232 352.38667 



15 



14 



13 



12 



568.32987 
492.55183 
422.67311 
358.64093 
300.39309 
11 



579.47499 513.05975 

502.26737 440.35862 

431.02140 373.63861 

365.69043 312.85154 

306.21917 257.93944 

252.54154 208.83194 
10 9 
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TABLE VII: Critical parameters f3 
V{r) = + 0.001 exp(r/ro)]. 



2mgrQh ^ of the Woods-Saxon potential 



n\l 







1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 



.050143914 .21236959 .45411555 

.40525072 .74663148 1.1653050 

1.0529125 1.5665339 2.1583382 

1.9863174 2.6709161 3.4345909 

3.2060640 4.0606118 4.9951932 

4.7125157 5.7361307 6.8408115 

6.5057952 7.6977523 8.9718549 

8.5859457 9.9456459 11.388595 

10.952985 12.479924 14.091224 

13.606920 15.300668 17.079885 

16.547758 18.407936 20.354690 

19.775499 21.801777 23.915725 

23.290146 25.482228 27.763062 

27.091699 29.449318 31.896761 13.822777 

31.180159 33.703074 12.256127 10.904356 

35.555527 10.683731 9.4306340 8.2552072 



.77426459 
1.6617334 
2.8284741 
4.2772382 
6.0094520 
8.0259370 
10.327215 
12.913644 
15.785486 
18.942937 
22.386154 
26.115262 
30.130362 



I. 1724439 1.6484948 
2.2359950 2.8880871 
3.5768855 4.4034793 
5.1986738 6.1987162 
7.1030599 8.2757388 
9.2910171 10.635665 

II. 763171 13.279221 
14.519959 16.206914 
17.561701 19.419123 
20.888645 22.916139 
24.500987 26.698194 
28.398884 16.962111 



2.2023382 
3.6179882 
5.3081667 
7.2772089 
9.5272608 
12.059576 
14.874975 
17.974033 
21.357183 
25.024758 



2.8339302 
4.4256753 
6.2908725 
8.4340224 
10.857441 
13.562505 
16.550125 
19.820942 
23.375440 



15.389823 13.842647 
12.372904 11.000819 
9.6302842 8.4339165 
7.1574531 6.1373746 



15 



14 



13 



12 



11 



18.543603 
15.318555 
12.373427 
9.7065447 
7.3152612 
5.1949761 
10 



20.137676 
16.804665 
13.753200 
10.982209 
8.4901060 
6.2743276 
4.3302634 
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FIG. 1: Function Fo(/3) = limr_).oo x'oi^) for the Yukawa potential {I = 0). Zeros of -Fo(/5) present 

critical parameters /3„. 
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FIG. 2: Function Fj{f3) = limr_j>oo '^('^) for the Gaussian potential (/ = 7). The abscissas of the 
points of intersection of t = Fj{f3) with hnes t = 6i — tt/2 + irn give the desired critical parameters 
/3 = /3„. 
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